Draft version February 2, 2008 

Preprint typeset using LATpjX style cmulateapj v. 25/04/01 



THE SZ EFFECT SIGNATURE OF EXCESS ENTROPY IN DISTANT, MASSIVE CLUSTERS 



m 

O 
O 



(N 
> 

in 

m 
O 
m 
o 

Of 

6 

C3 



13 



Ian G. McCarthy 1,2 , Gilbert P. Holder 3 , Arif Babul 1,4 , and Michael L. Balogh 5 

department of Physics & Astronomy, University of Victoria, Victoria, BC, V8P 1A1, Canada 
3 School of Natural Sciences, Institute for Advanced Study, Princeton, NJ, 08540, USA 
5 Department of Physics, University of Durham, South Road, Durham, DH1 3LE, UK 
Draft version February 2, 2008 

ABSTRACT 

Studies of cluster X-ray scaling relations have led to suggestions that non-gravitational processes, e.g., 
radiative cooling and/or "preheating", have significantly modified the entropy of the intracluster medium 
(ICM). For the first time, we test this hypothesis through a comparison of predicted thermal Sunyaev- 
Zeldovich (SZ) effect scaling relations with available data from the literature. One of the relations that 
we explore, in principle, depends solely on SZ effect observations, thus offering an X-ray independent 
probe of the ICM. A detailed comparison of the theoretical relations with the largest compilation of 
high-z SZ effect data to date indicates that the presence of an entropy floor is favored by the data. 
Furthermore, the inferred level of that floor, K$ > 300 keV cm 2 , is comparable to that found in studies 
of X-ray scaling relations of nearby massive clusters. Thus, we find no evidence for significant evolution 
of the entropy floor out to z ~ 0.7. We further demonstrate that the high quality data to be obtained 
from the upcoming Sunyaev-Zeldovich Array (SZA) and the (soon-to-be) upgraded Owens Valley Radio 
Observatory (OVRO) array will open powerful new windows into the properties of the ICM. Specifically, 
the new measurements will allow for accurate measurements of the ICM entropy for even the most distant 
galaxy clusters. 

Subject headings: cosmology: cosmic microwave background — cosmology: theory — galaxies: clusters: 
general — X-rays: galaxies: clusters 



1. INTRODUCTION 

The failure of theoretical self-similar X-ray scaling re- 
lations to match observed trends has led to suggestions 
that important non-gravitational processes, such as ra- 
diative cooling and/or "preheating", are significantly af- 
fecting the structure and appearance of the intracluster 
medium (ICM). Models and simulations of clusters which 
attempt to explicitly take into account the effects of cool- 
ing and/or heating produce clusters which have higher 
mean entropies ("excess" entropy) than those produced 
by models which neglect these processes. In some cases, 
the model clusters possess cores in their entropy profiles 
commonly referred to as the "entropy floor" . The presence 
of this entropy floor, in turn, modifies the predicted X-ray 
scaling relations, bringing them in to much closer agree- 
ment with the observed correlations (e.g., Kaiser 1991; 
Evrard & Henry 1991; Bower 1997; Balogh et al. 1999; Wu 
et al. 2000; Bryan 2000; Tozzi & Norman 2001; Borgani et 
al. 2001; Voit & Bryan 2001; Babul et al. 2002; McCarthy 
et al. 2002; Thomas et al. 2002; Voit et al. 2002; Dave et 
al. 2002; Lloyd-Davies et al. 2003). Direct observational 
evidence for an "entropy floor" in nearby groups and low 
mass clusters has been presented by Ponman et al. (1999) 
and Lloyd-Davies et al. (2000). 

To date, only X-ray observations have yielded informa- 
tion about the entropy floor but because the X-ray surface 
brightnesses of groups/clusters suffer the effects of cosmo- 
logical dimming [the bolometric surface brightness scales 
as (l+z)~ 4 ] , both direct and indirect studies of the entropy 
floor via X-ray observations have generally been limited to 
low redshift (z < 0.2) systems. An independent test of the 
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entropy floor hypothesis which could also provide informa- 
tion on high redshift clusters and, therefore, the evolution 
of the non-gravitational processes that give rise to the en- 
tropy floor, would be extremely useful. In a companion 
paper (McCarthy et al. 2003, hereafter MBHB03), we ar- 
gued that a number of scaling relations based entirely or in 
part on thermal Sunyaev-Zeldovich effect (Sunyaev & Zel- 
dovich 1972; 1980; hereafter, referred to as the SZ effect) 
observables can, potentially, be used for both of these pur- 
poses. Our analysis indicated that even current SZ effect 
observational data from, for example, the Berkeley Mary- 
land Illinois Association (BIMA) and Owens Valley Radio 
Observatory (OVRO) arrays and the Ryle Telescope, when 
compared to the predicted correlations, can be used to tell 
us something about the entropy floors of distant, massive 
clusters. 

The primary focus of the present paper is a compari- 
son of these predicted scaling relations with available data 
from the literature to determine if the SZ effect data fa- 
vor the existence of an entropy floor and how the inferred 
level of the entropy floor compares with that required to 
explain local X-ray trends (which require Kq > 300 keV 
cm 2 for massive clusters; Tozzi & Norman 2001; Babul et 
al. 2002; McCarthy et al. 2002). This is the first time that 
the SZ effect has been used as a probe of non- gravitational 
entropy in galaxy clusters. 

A number of new SZ experiments, which will greatly 
improve the quality of the observations, are being planned 
or are already under construction. An additional goal of 
the present study, therefore, is to examine the efficacy of 
two of these experiments, the Sunyaev-Zeldovich Array 
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(SZA) and the (soon to be) upgraded OVRO array, to con- 
strain the properties of the excess entropy in distant clus- 
ters. By constructing and analysing "mock" observations 
which take explicit account of the expected instrumental 
responses of these arrays, we quantify the accuracy with 
which the level of the entropy floors of distant clusters can 
be inferred by future data as a function of rcdshift. 

The present paper is outlined as follows: in §2, we dis- 
cuss and select available data from the literature; in §3, 
we compare these data to our theoretical scaling relations; 
in §4, we assess the ability of the upcoming SZA and up- 
graded OVRO array to probe the entropy floors of distant 
clusters; and in §5; we discuss and summarize our results. 
We assume VL m = 0.3, £1a = 0.7, and H = 75 km s _1 
Mpc -1 and work in physical units (e.g., M Q rather than 
Ii~ 1 Mq) throughout the paper. 

2. OBSERVATIONAL DATA 

Observations of the SZ effect have advanced tremen- 
dously over the last decade or so and routine high signal- 
to-noise measurements of the effect are now being made 
with a variety of instruments (including both single-dish 
and interferometric experiments) at a variety of wave- 
lengths (spanning radio down to the submillimeter) . There 
now exist published data for some 30-40 clusters (e.g., 
Jones et al. 1993; Herbig et al. 1995; Carlstrom et al. 
1996; Myers et al. 1997; Holzapfel et al. 1997; Hughes & 
Birkinshaw 1998; Pointecouteau et al. 1999; 2001; 2002; 
Patel et al. 2000; Holzapfel et al. 2000; Joy et al. 2001; 
Mason et al. 2001; Grego et al. 2001; Grainge et al. 2002a; 
Grainge et al. 2002b; Jones et al. 2003; Reese et al. 2002; 
Cotter et al. 2002; LaRoque et al. 2003; Cantalupo et al. 
2003). With this influx of new data, the sample is large 
enough to make statistically significant comparisons be- 
tween observed and theoretically predicted SZ effect scal- 
ing relations. Such comparisons test our understanding 
of the ICM and clusters in general. In this section, we 
compile and discuss cluster SZ effect data from the liter- 
ature with the intention of comparing it to theoretically 
predicted scaling relations in §3. 

Of the 30-40 clusters that have published SZ effect data, 
we are particularly interested in those clusters which lie at 
high redshift. As already mentioned, the vast majority of 
studies on the X-ray scaling relations of clusters have been 
for low rcdshift (z < 0.2) systems but information on the 
entropy floors of more distant clusters is scant. Because 
the SZ effect is not subject to cosmological dimming, SZ ef- 
fect scaling relations potentially offer a way of probing the 
non-gravitational entropy of even the most distant galaxy 
clusters. 

A search of the literature for high redshift clusters with 
SZ effect data yields 22 different clusters in the range 
0.14 < z < 0.78, many of which were observed multiple 
times. The clusters are listed in Table 1 along with their 
rcdshifts, central Compton parameters (yo), frequency- 
independent integrated SZ effect flux densities within the 
central 1' (S v , arc / '/„), and the bibliographic references (ref- 
erences in parentheses indicate that the cluster was ob- 
served multiple times). We are particularly interested in 
yo and S v ^ arc / '/„ because a number of scaling relations 
based on these two quantities are expected to be quite 
sensitive to the entropy floor level of galaxy clusters (see 



MBHB03). 

Below, we discuss how we extract yo and S u . arc / '/„ from 
the observational data. At present, it is not possible to 
directly measure either of these quantities. Observations 
of the SZ effect filter large-scale emission while finite res- 
olution smears out small-scale structures. As discussed in 
MBHB03, fitting a model (such as the isothermal (3 model; 
Cavaliere & Fusco-Femiano 1976; 1978) to the SZ effect 
data provides a method for effectively removing these ef- 
fects and estimating yo and S v . arc j f v , but it should be 
kept in mind that such quantities are inferred and model- 
dependent. Provided the smallest angular scale resolved 
is comparable to the typical scale over which the clus- 
ter varies, the estimated central Compton parameter will 
be reliable, while inferred flux densities will be suspect 
when extrapolated beyond the filtering scale of the ob- 
servations. For current interferometric observations [such 
as those obtained with the BIMA and OVRO arrays and 
the Ryle telescope], the highest angular resolution for SZ 
measurements is typically smaller than the core radius of 
the cluster (~ 30") while the large-scale filtering normally 
becomes important on scales larger than about 2'. There- 
fore, it can be expected that the inferred values of yo and 
Sv,arc/ fv (the flux density within the central 1') should be 
accurate. 

To calculate the total frequency-independent SZ effect 
flux density within the central 1' for each of these clusters, 
we must reconstruct each of the SZ effect "surface bright- 
ness" profiles, y{9). The majority of the clusters listed 
in Table 1 were modeled using the spherical isothermal (3 
model. In this model, the ICM is assumed to be isothermal 
and has a density distribution described by 

/ r 2\ -3/5/2 
n e {r) = n e o[l + (1) 

were n e o is the central electron density, r c is the cluster 
core radius, and (3 is the power-law index. This leads to a 
SZ effect surface brightness profile 

/ f)2\ 1/2-3/3/2 

vV)=Vo(l + ¥ J (2) 

where y(9) is the Compton parameter evaluated at a pro- 
jected position 9 = r/D a from the cluster center and is 
proportional to the integrated pressure along the line-of- 
sight through the cluster. In addition, 9 C = r c /D a (D a is 
the angular diameter distance). 

The central Compton parameters, yo, are converted 
from the central SZ effect temperature decrements (ref- 
erences are given in col. 5) using the relation 

AT o fx \ 

where ATo is the central SZ effect temperature decrement 
and x = hv/kTcMB is the dimensionless frequency (Tcmb 
is the temperature of the present-day cosmic microwave 
background — 2.728 K; Fixscn et al. 1996, and v is the 
observing frequency). We ignore the complication of rela- 
tivistic effects, which only modify the Compton parameter 
of the hottest clusters by a few percent (Itoh et al. 1998; 
Nozawa et al. 2000). 
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We use the best-fit [3 model parameters [yo (ATb), /?, 
and 9 C ] from the literature (discussed below) to reconstruct 
y{9). The surface brightnesses are then numerically inte- 
grated within the central I 1 (the result is symbolized by 
Vint) via 



0—1 

y in t(< e = 1') = 2tt / y(9')fl'^' 
Jo 



(4) 



and, finally, converted into a frequency-independent flux 
density through 



Sv^arc/ fv — y% 



2(kT ( 



CMB) 



(h c y 



(5) 



where f v is a function of the dimensionless frequency x 
and is defined in MBHB03. 

Ideally, the three (3 model parameters would be deter- 
mined by fitting the model to the SZ effect data. How- 
ever, current SZ effect data cannot tightly constrain these 
parameters when all three are left to vary independently 
(e.g., Carlstrom ct al. 1996; Grego et al. 2000; 2001; 
Pointecouteau et al. 2001; 2002). This problem is often 
circumvented by adopting the best-fit values of (3 and 9 C 
(the "shape" parameters) determined from fitting to the 
X-ray surface brightness profile of the cluster and leaving 
only the normalization, yo, to be determined from fitting 
to the SZ effect data (e.g., Pointecouteau et al. 1999; 2001; 
2002; Patel et al. 2000; Jones et al. 2003; LaRoque et al. 
2003; Grainge et al. 2002b). In the case of clusters with 
moderate redshifts (which is a good description of most 
of the clusters in Table 1), X-ray data still provide better 
constraints on (3 and 9 C than do SZ effect data. Thus, 
in the cases of A1914 and RXJ2228+2037, we use the 
X-ray-determined values of these parameters to calculate 
S v ,arc/ fv However, a better approach is to use all avail- 
able data (both the SZ effect and X-ray data) to constrain 
these parameters. In estimating Hubble's constant from a 
sample of 18 distant clusters, Rccsc ct al. (2002) did just 
this. The shape parameters, central X-ray surface bright- 
ness, and yo were all determined simultaneously by using 
a joint maximum-likelihood analysis of both SZ effect and 
X-ray data. Because their sample is large (in fact, it is 
the largest sample of SZ effect clusters observed to date) 
and homogeneously analysed and their method takes ad- 
vantage of all ICM imaging data, we preferentially use the 
values of yo, (3, and 9 C measured by Reese et al. (2002) 
when multiple measurements for a particular cluster are 
available (which is the case for roughly half of the clusters 
listed in Table 1). We note that the agreement between 
different studies is reasonably good (estimates of yo differ 
by < 20% from study to study, e.g., Holzapfel et al. 1997; 
Jones et al. 2003). Using the best-fit parameters of Reese 
et al. (2002), we calculate S v ^rcl 1 fv for these 18 clusters. 
Published values for the shape parameters of the two re- 
maining clusters in Table 1, A2204 and Zwicky 3146, are 
not available. 

To calculate the uncertainty associated with the SZ ef- 
fect flux density, we vary yo, (3, and 9 C within their al- 
lowable ranges (each of the three free parameters has an 
associated statistical uncertainty) to determine the max- 
imum and minimum possible flux density of the cluster. 
This method actually overestimates the uncertainty asso- 



ciated with the flux density as there is a known correlation 
between the (3 and 9 C parameters for current SZ effect data 
(e.g., Grego et al. 2000; 2001). 

We are also interested in mapping out correlations be- 
tween SZ effect and X-ray observables. In Table 2, we 
list the total dark matter masses within r$oo [M(r§oo)]i 
the mean emission- weighted gas temperatures (Tx), and 
bolometric X-ray luminosities (Lx) of these clusters. We 
also list whether the cluster has a sharp centrally-peaked 
X-ray surface brightness profile which presumably indi- 
cates of the presence of a cooling flow (CF indicates the 
cluster is a "cooling flow cluster" , NCF indicates that it is 
not; Allen 2000; Reese et al. 2002) and the references for 
M(r5oo), Tx, and Lx (respectively). Below, we discuss 
the X-ray data in columns 2-4 of Table 2. 

For the total dark matter masses within r 50 o (col. 2), we 
turn to the study of Ettori & Fabian (1999). By construct- 
ing X-ray surface brightness profiles that are based on the 
"universal" dark matter density (NFW) profile (Navarro 
et al. 1997) and comparing these to data from ROSAT, 
these authors were able to model the underlying dark mat- 
ter density for nine of the clusters listed in our Table 1. 
We use their best-fit parameters (listed in their Tables 1 
and 2) to recover the dark matter density profiles of these 
clusters. We take care to properly scale these parameters 
for our assumed cosmology. We then integrate the density 
profiles out to r 50 o to determine the total dark matter mass 
within that radius. We are unable to estimate the uncer- 
tainty on the masses, as there is no reported uncertainty 
for the best-fit NFW parameters of these clusters. 

The most recently determined values for the ICM tem- 
peratures of these clusters are also listed in Table 2 (col. 
3). Of the 21 clusters with reported temperatures, 16 
have temperatures determined via fitting to ASCA X-ray 
spectral data (Allen 2000; White 2000; Novicki, Sornig, & 
Henry 2002; Jones et al. 2003) and five have measurements 
based on fits to new Chandra X-ray spectral data (Marke- 
vitch & Vikhlinin 2001; Machacek et al. 2002; Vikhlinin 
et al. 2002). The cluster RXJ2228+2037 docs not have 
a temperature deduced from X-ray spectral analysis, al- 
though its central temperature was estimated by a com- 
bined analysis of X-ray and SZ effect imaging data (Pointe- 
couteau et al. 2002). In the interest of homogeneity, how- 
ever, we do not include this cluster in our analysis of scal- 
ing relations involving the emission-weighted gas tempera- 
tures because it is unclear how the reported temperature is 
related to the X-ray emission-weighted temperature. For 
the clusters that have sharp centrally-peaked X-ray surface 
brightnesses and apparently harbor massive cooling flows 
(Allen 2000; Reese et al. 2002; see col. 5), we use cool- 
ing flow-corrected temperatures. The temperatures were 
corrected by fitting the X-ray spectra with a multi-phase 
plasma model that explicitly takes into account the cooler 
emission from the cluster core (Allen 2000). We note that 
for the clusters in this sample, cooling flow correction only 
slightly increases the temperature of a cluster (by about 
10% or 1 keV) and does not significantly affect our results. 

Finally, column 4 lists the bolometric X-ray luminosi- 
ties of these clusters. We preferentially select data that 
has published uncertainties (e.g., White, Jones, & Forman 
1997; Novicki, Sornig, & Henry 2002). For data without 
published uncertainties, we select on the basis of X-ray 
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satellite in the following order: Chandra/ 'XMM-Newton, 
ASCA, and RO SAT /Einstein. Of the 17 clusters that 
have published bolometric luminosities, seven were deter- 
mined using Einstein imaging data (White, Jones, & For- 
man 1997), one was determined using ROSAT imaging 
data (Allen & Fabian 1998), seven were determined using 
ASCA imaging data (Mushotzky & Scharf 1997; Novicki, 
Sornig, & Henry 2002), and two were determined using 
new Chandra imaging data (Vikhlinin et al. 2002). Unfor- 
tunately, cooling flow-corrected luminosities (which nor- 
mally entails excising the central 100-200 kpc of the CF 
cluster; e.g., Markevitch 1998; Vikhlinin et al. 2002) were 
not available for the CF clusters in Table 2. 

Tables 1 and 2 comprise the largest compiled sample of 
high rcdshift SZ effect clusters to date. We note that the 
current compilation is not statistically complete or homo- 
geneous, as discussed in Reese et al (2002). Clusters with 
bright radio point sources were avoided, for the most part, 
and targets were selected primarily on the basis of X-ray 
luminosity. This compilation can not be simply defined 
as being flux-limited or luminosity-limited, for example, 
and it is not complete over the range of redshifts studied. 
Without a simple selection function it is difficult to assess 
possible systematic effects that could be introduced by not 
having a sample that adequately reflects the demographics 
of the general population of galaxy clusters. 

Below, we compare the data in Tables 1 and 2 to our 
theoretically predicted scaling relations (see Table 1 of 
MBHB03) and attempt to determine if the data require 
an entropy floor (similar to that of X-ray data of nearby 
clusters) and, if so, what is the level of that entropy floor. 

3. RESULTS 

3.1. Comparing Theory to Observations 

A detailed analysis of how the presence of excess en- 
tropy alters SZ effect scaling relations was presented in the 
companion paper MBHB03. These relations were derived 
using the entropy injection (preheated) models developed 
by Babul et al. (2002). Following standard practice, we 
fitted simple power-law models to the theoretical scaling 
relations. Since a number of the relations were quite sen- 
sitive to rcdshift and existing SZ effect data span a wide 
range in z (see Table 1), we fitted power-law models which 
were an explicit function of both redshift and entropy floor 
level (K ) to the theoretical relations; i.e., for arbitrary 
cluster parameters X and Y (e.g., y and Tx) 

log 10 y = a(K ,z)log 10 X + b(K Q ,z) (6) 

With these relations, it is possible to quickly and accu- 
rately compute how various cluster properties scale with 
the SZ effect of the Babul et al. (2002) model clusters at 
any rcdshift between 0.1 < z < 1 and with any entropy 
floor in the range 100 keV cm 2 < Kq < 700 keV cm 2 (see 
the Table 1 of MBHB03). This is the only way a fair com- 
parison can be made between the existing data and the 
theoretical models. 

We compare the theoretical scaling relations to the ob- 
servational data compiled in Tables 1 and 2 with a x 2 
statistic 



2 = v [lQSYi-a(K Q ,z)logXi-b(K Q ,z)] 2 
X h * 2 Yi +a(K ,zro- Xi [l) 

where u\. and a\. are the measurement errors in log-space 
and are estimated by multiplying the measurement errors 
in linear-space by the weighting factors (liloge) -1 and 
pQloge) -1 , respectively (see, e.g., Press et al. 1992). 

We specify the optically-determined redshifts for each 
of the clusters and then determine the best-fit value of 
entropy floor level, Kq, by minimizing the x 2 (assuming 
that the value of Kq is the same for all clusters). Unless 
stated otherwise, quoted error bars are for the 95.4% (2a) 
confidence level (corresponding to A% 2 = 4). We assess 
the quality of the fit by calculating the reduced-x 2 (x 2 )- 
For the self-similar model (whose scaling relations depend 
on z only; i.e., there arc no free parameters), we calculate 
the x 2 an d xl and compare this with that of the best-fit 
entropy floor model. 

Below, we examine seven different scaling relations in- 
volving the two primary SZ effect observables (the cen- 
tral Compton parameter and the integrated flux density). 
First, we explore the trend that exists between the two SZ 
effect observables. This is an extremely interesting test 
since it can potentially be measured entirely independent 
of X-ray observations. In fact, in §4, we show that future 
data from the SZA and the upgraded OVRO array will 
allow one to constrain the central entropy distribution of 
clusters all the way out to z ~ 2 using this relation. The 
next two trends that we study, between the two SZ effect 
observables and the mass of a cluster, are also interesting 
because they too can potentially be measured independent 
of X-ray observations. Aside from X-ray observations, one 
can estimate the mass of a cluster via gravitational lensing 
or through galaxy velocity dispersion measurements. We 
expect that the yo — M(rsoo) relation, in particular, will 
be very sensitive to the presence of excess entropy. Fi- 
nally, the four remaining correlations that we look at are 
between the two SZ effect observables and the two primary 
X-ray observables, i.e., the mean emission- weighted tem- 
perature (Tx) and the bolometric luminosity (Lx)- The 
correlations between the SZ effect observables and the clus- 
ter temperature are expected to be excellent probes of the 
central entropy of clusters and, as such, we discuss these 
ahead of those involving Lx ■ 

3.2. A SZ Effect Only Relation 

We start by examining the correlation between the two 
SZ effect quantities, S u ^ ar c/ fv and y (see Fig. 2 of 
MBHB03). Even though theoretical arguments suggest 
that this relation is not as sensitive to the presence of 
excess entropy as some of the other relations we will dis- 
cuss later, it is by far the most interesting. As we already 
mentioned, this relation can potentially be measured en- 
tirely through SZ effect observations, offering a completely 
X-ray-independent probe of the intracluster gas. This is 
discussed in more detail below. 

Fitting all 20 clusters in Table 1 for which we have esti- 
mates of both S Ut arc/ fv and yo, we find a best-fit entropy 
floor level of K = 540t\l° keV cm 2 with a xl = 37.77/19 
= 1.99. A residual plot demonstrating the quality of our 
best-fit entropy floor model is presented in Figure 1 (top 
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Fig. 1. Residual plots for the Su^arc/fu — yo relation. Top: 
Residuals from a comparison between the entropy floor model with 
K = 540 kcV cm 2 and the data in Table 1 (x 2 = 1-99). Mid- 
dle: Neglecting A1835 and RXJ1347.5-1145 - for reasons discussed 
in the text - greatly improves the fit: Kq = 570 keV cm 2 with 
X 2 = 0.99. Bottom: Residuals from a comparison between the self- 
similar model and the data in Table 1 (x 2 = 12.21). Filled symbols 
represent low redshift (z < 0.3) clusters, open symbols represent high 
rcdshift (z > 0.3) clusters, triangles represent cooling flow clusters, 
and squares represent non cooling flow clusters. For clarity, error 
bars for the abscissa are not displayed. 

panel). The residuals from a comparison between the self- 
similar model and the data in Table 1 is presented in the 
bottom panel {xl = 244.15/20 = 12.21). Note that the 
residuals in the top panel are generally consistent with 
the zero line (indicating a relatively good fit), whereas 
the residuals in the bottom panel are systematically too 
high, indicating that the self-similar model predicts val- 
ues of S v . arc l f v (for a fixed yo) that are systematically 
lower than observed. This is equivalent to saying that 
the self-similar model predicts clusters that have SZ effect 
surface brightness profiles that are too centrally-peaked. 
The entropy floor models are able to match the observa- 
tions much better because the addition of an entropy floor 
flattens the pressure profiles of clusters (e.g., Fig. 1. of 
MBHB03) and since the SZ effect is proportional to in- 
tegrated line-of-sight pressure of the cluster, this leads to 
flatter surface brightness profiles. 

There is little doubt that the best-fit entropy floor model 
provides a much better fit than the self-similar model to 
the observational data, but the best-fit entropy floor model 
itself does not provide a statistically acceptable fit to the 
data (note the high value of the xl)- However, the resid- 
uals in the top panel of Figure 1 clearly show two obvious 
outliers, A1835 and RXJ1347.5-1145. Neglecting these two 
clusters, the quality of the fit to the whole sample is sig- 
nificantly improved: we find Kq = 575^^5 keV cm 2 with 
xl = 16.79/17 = 0.99. The residuals for this fit are plotted 
in the middle panel of Figure 1. This result is consistent 



with that determined from studies of X-ray scaling rela- 
tions of nearby massive clusters, which require Kq > 300 
keV cm 2 (e.g., Tozzi & Norman 2001; Babul et al. 2002; 
McCarthy et al. 2002). 

It is interesting that the two outliers, A1835 and 
RXJ1347. 5-1145, have extraordinarily large cooling flow 
mass deposition rates. Both apparently deposit several 
thousand solar masses of gas each year and are among the 
most massive cooling flow clusters known (Allen 2000). In 
fact, a closer inspection of the residuals (particularly in 
the left-hand panel) reveals that there is a small system- 
atic difference between cooling flow clusters (triangles) and 
non cooling flow clusters (squares), even if one neglects 
these two outliers. While in principle the SZ effect should 
be less susceptible than the X-ray emission to the effects 
of cooling flows (because SZ effect is proportional to n e , 
whereas the X-ray emission is proportional to n 2 ), recall 
that both of the SZ effect quantities used in Figure 1 were 
inferred through fitting to both X-ray and SZ effect data. 
Thus, the strong central surface brightness peaks present 
in the X-ray images of these cooling flow clusters will have 
an impact on the implied SZ effect observables. Excluding 
A1835 and RXJ1347.5-1145, however, there is no statis- 
tical evidence for a difference in Kq when fitting to the 
cooling flow and non cooling flow clusters separately or 
to the whole sample. Therefore, it appears that only the 
most extreme cooling flow clusters could have significantly 
different entropy histories. 

Because the sample in Table 1 is reasonably large and 
spans a wide range of redshifts, it is possible to use 
the catalog to get some idea of how the entropy floor 
level of clusters evolved with cosmic time. We split the 
sample up into two large redshift bins: (1) "low" red- 
shift clusters (z < 0.3), and (2) "high" rcdshift clusters 
(z > 0.3). Fitting only the low redshift clusters, we derive 
K = 5051^0 keV cm2 (xl = 10.46/10 = 1.05), which is 
consistent with the value derived from fitting to the en- 
tire sample. The high redshift clusters, however, prefer a 
slightly elevated entropy floor with Kq — 640^215 keV cm 2 
(xl = 5.84/6 = 0.97), although the difference between the 
low and high redshift clusters is not statistically signifi- 
cant. Thus, there is no good evidence that K evolves 
significantly out to z <~ 0.7, at least on the basis of this 
test. A stronger test of this hypothesis will soon be pos- 
sible, as the list of high redshift clusters observed through 
the SZ effect is growing rapidly (Carlstrom and Joy and 
collaborators, for example, have now made detections in 
21 clusters with z > 0.45; Reese et al. 2002). 

Finally, to what extent the present S u ^ arc / f u — yo rela- 
tion studied here is independent of previous X-ray results 
is debatable. Both S v ^ arc / /„ and yo were inferred by fit- 
ting a model to the SZ effect surface brightness profiles 
of clusters and, furthermore, this surface brightness model 
(the isothermal (3 model) had two of its three free param- 
eters constrained to be the same as that deduced from 
X-ray observations (Jones et al. 2003; Pointecouteau et 
al. 2002) or from X-ray and SZ effect observations (Reese 
et al. 2002). Yet, it is also clear that this test is different 
from any other scaling relation examined to date. Ideally, 
both S v . arc l '/„ and yo could be measured directly or, fail- 
ing that, determined from fitting a model to only the SZ 
effect data. Current SZ effect data alone, however, cannot 
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Fig. 2. The observed and predicted j/o — M(rsoo) relations. The 
symbols have the same meaning as in Figure 1. The dotted, short- 
dashed, long-dashed, dot-dashed, and solid lines represent the self- 
similar and Kq = 100, 300, 500, and 700 keV cm 2 entropy floor 
models, respectively. 

tightly constrain these quantities. In §4, we show that the 
high quality data to be produced by the upcoming SZA 
and the (soon to be) upgraded OVRO array will make it 
possible to accurately measure these quantities for mas- 
sive clusters at virtually any redshift and without having 
to use any X-ray results. 

3.3. The SZ effect-M(r 500 ) Relations 

The next set of scaling relations we examine are the SZ 
effect - cluster mass relations. Theoretical arguments sug- 
gest that these relations should be very sensitive to the 
presence of entropy floor, at least when the SZ effect is 
measured near the cluster center (e.g., yo)- Interestingly 
these trends too can also potentially be measured indepen- 
dent of X-ray results. Future SZ effect observations will 
allow one to estimate y and S„. arc /f u accurately purely 
through SZ effect surface brightness profiles (see §4), while 
both strong and weak lensing are increasingly being used 
to measure the mass profiles of clusters out to radii com- 
parable in size to that of r 50 o (e.g., Clowe & Schneider 
2001). At present, however, only a few of the clusters in 
our sample have been weighed using lensing. 

In Figure 2, we plot the observed yo — M(rsoo) rela- 
tion. This is superimposed on the predicted z — 0.2 
yo — M(r 500 ) relations for the self-similar model (dotted 
line) and the Kq = 100 (short-dashed), 300 (long-dashed), 
500 (dot-dashed), and 700 keV cm 2 (solid) entropy floor 
models. Because the majority of the clusters with pub- 
lished masses in Tables 1 and 2 lie in a narrow redshift 
range around z ~ 0.2, it is possible to qualitatively com- 
pare the theoretical models to the data "by eye" before 
using the quantitative method outlined in §3.1. A simple, 
neat, and fair qualitative comparison was not possible for 
the previous relation (S Vtarc / f v — yo), as the data spanned 
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Fig. 3. Residual plots for the J/o — -^(T500 ) relation. Left: Residuals 
from a comparison between the entropy floor model with Ko = 500 
keV cm 2 and the observational data (x 2 = 1.00). Right: Residuals 
from a comparison between the self-similar model and the observa- 
tional data (x 2 = 41.80). The symbols have the same meaning as in 
Figure 1. For clarity, error bars for the abscissa are not displayed. 

a wide range of redshifts and because that relation is es- 
pecially sensitive to z (S u scales as 1/-D 2 ). 

By visual inspection of Figure 2, it is obvious that only 
the high entropy floor models (Ko > 300 keV cm 2 ) pro- 
vide a reasonable fit to the observational data. In addition, 
the observed correlation does not seem to depend on the 
cooling flow status or redshift of the clusters, although the 
sample is far too small to make any robust conclusions to 
this effect. Fitting all nine clusters with the method out- 
lined in §3.1, our best-fit entropy floor level is K a = 500+gg 
keV cm 2 with xl = 8.02/8 = 1.00. This is consistent with 
the results derived in §3.2 and with X-ray observations of 
nearby massive clusters. A plot of the residuals between 
the data and the K — 500 keV cm 2 model is shown in the 
left-hand panel of Figure 3. Also shown (right-hand panel) 
are the residuals of a comparison between the data and the 
self-similar model (xl = 376.05/9 = 41.80). The residuals 
for the Ko = 500 keV cm 2 model display a tight scatter 
about the zero line, while the residuals for the self-similar 
model indicate that yo is observed to be much lower [for a 
fixed value of M(rsoo)] than predicted by this model. The 
entropy floor models with K ~ 500 keV cm 2 are able to 
provide a good match to the data because the addition of 
an entropy floor reduces the gas pressure near the centers 
of clusters (see MBHB03 for a detailed discussion). This, 
in turn, reduces the magnitude of y a - The mass within 
rsoo, however, is unaffected by the modification of the gas 
entropy. 

Although the available yo — M(r 500) data exhibit only 
a very small amount of scatter about the K ~ 500 keV 
cm 2 relation (and the xl indicates a very good fit), the 
estimated error bars on our best-fit value of Ko from this 
relation are almost certainly too small. We say this be- 
cause (1) we were unable to calculate any uncertainty for 
-^(^soo) as there were no published error bars for the best- 
fit NFW parameters for the clusters in Figures 2 and 3, 
and (2) the sample is too small to get any kind of a han- 
dle on the systematic errors associated with the cluster 
masses. For example, estimates of cluster masses from 
gravitational lensing very often differ from those deter- 
mined from X-ray data (sometimes by up to a factor of 
two; e.g., Miralda-Escude & Babul 1995; Wu & Fang 1997; 
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Fig. 4. The observed and predicted yo — Tx relations. The symbols 
have the same meaning as in Figure 1. The dotted, short-dashed, 
long-dashed, dot-dashed, and solid lines represent the self-similar 
and K = 100, 300, 500, and 700 keV cm 2 entropy floor models, re- 
spectively. Thin lines arc for z = 0.2 and thick lines are for z = 0.5. 
As discussed in the text, the models with Kg ~ 300 keV cm 2 provide 
the best description of the observations. 

Ota et al. 2002) and it is not yet clear why this happens 
to be the case. There could be a systematic problem with 
the X-ray determined masses (e.g., Allen 1998). 

Aside from the yo — M(rsoo) relation, we also explore 
the 

Su.arc/ fu ~ M^r^oo) trend. This trend could provide 
a consistency check of the results discussed immediately 
above. Unfortunately, the available SZ effect flux density 
- mass data do not constrain K$. This is not completely 
unexpected since the S v . arc / f v — M(r§oo) relation is much 
less sensitive than the yo — M(r 50 o) relation to the en- 
tropy floor level of clusters (MBHB03). Future data from, 
for example, the SZA will allow for much more precise de- 
terminations of the SZ effect flux densities of clusters and 
we expect that the future S v . arc j f v — M(rsoo) relation will 
place much tighter constraints on K - 

3.4. The SZ effect— Tx Relations 

Plotted in Figure 4 is the observed yo — Tx relation. In 
addition to the observational data, we have also plotted 
the z = 0.2 (thick lines) and z = 0.5 (thin lines) theoret- 
ical relations. The various line types (e.g., solid, dotted) 
have the same meanings as in Figure 2. Again, because the 
theoretical relations are not a strong function of redshift, 
a preliminary visual comparison of the data and models is 
feasible. 

First, it is clear that the data indicate that there is a 
fairly tight correlation between a cluster's central Comp- 
ton parameter and its emission-weighted gas temperature. 
This correlation seems to hold true irrespective of whether 
the clusters have cooling flows or not. This may be ex- 
pected since we have used cooling flow-corrected tempera- 
tures. However, even without this correction the data still 




T x (keV) 

Fig. 5. Residual plots for the yo — Tx relation. Left: Residuals 
from a comparison between the entropy floor model with Ko = 300 
keV cm 2 and the observational data (x 2 = 1.32). Neglecting A1835 
and RXJ1347.5-1145, we find a best fit of K = 305 keV cm 2 with xl 
= 1.24. Right: Residuals from a comparison between the self-similar 
model and the observational data (x 2 = 21.41). The symbols have 
the same meaning as in Figure 1. For clarity, error bars for the 
abscissa arc not displayed. 

exhibit a tight trend. For the clusters in this sample, cool- 
ing flow correction only slightly increases the temperature 
of a cluster (by about 10% or 1 keV) and, therefore, does 
not significantly affect the results. There is also a hint of 
a slight systematic difference in the yo — Tx relations for 
low and high redshift clusters (filled and open symbols, 
respectively). In particular, the relation for the high red- 
shift clusters has a normalization that is slightly higher 
than that of the relation for low redshift clusters (ignoring 
RXJ1347, whose temperature is quite uncertain). This is 
discussed in more detail below. 

Comparing the various theoretical relations to the obser- 
vational data, it is apparent that the standard self-similar 
model predicts a yo — Tx relation which is a poor match 
to the data. In particular, the normalizations of the self- 
similar relations arc roughly 2.5 times larger than observed 
(over the range 6 keV < Tx < 12 keV). The entropy floor 
models with K ~ 300 keV cm 2 , on the other hand, pro- 
vide a very good qualitative fit to the data. In MBHB03, 
we found that injecting the ICM with non-gravitational en- 
tropy tends to decrease yo and, at the same time, increase 
Tx- This reduces the normalization of the predicted trend 
between these two quantities and, as is apparent from Fig- 
ure 4, brings close agreement between theory and obser- 
vations. 

Fitting all 21 clusters in Tables 1 and 2 via the method 
outlined in §3.1, we obtain a best-fit entropy floor of 
Ko = 300t^ keV cm 2 with xl = 26.36/20 = 1.32. The 
residuals of a comparison between the K = 300 keV cm 2 
entropy floor model and the self-similar model with the 
observation data arc plotted in Figure 5. Note the strik- 
ing separation between the residuals of the two models. 
Theoretically, the yo — Tx relation is the most sensitive to 
the entropy floor level of any and all scaling relations that 
we have explored. 

If we exclude A1835 and RXJ1347.5-1145 from the fit, 
the entropy floor is essentially unchanged (K a = 305+g" 
keV cm 2 ) but the fit is improved (xl = 22.37/18 = 1.24). 
In addition, we verify that splitting the sample by cooling 
flow status does not significantly modify the best-fit value 
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Fig. 6. The observed and predicted yo — Lx relations. The symbols 
have the same meaning as in Figure 1. The dotted, short-dashed, 
long-dashed, dot-dashed, and solid lines represent the self-similar 
and Kq = 100, 300, 500, and 700 keV cm 2 entropy floor models, re- 
spectively. Thin lines are for z = 0.2 and thick lines are for z = 0.5. 

of Kq. As mentioned above, there appears to be a slight 
difference in the normalizations of the yo — Tx relations 
for low and high rcdshift clusters. Is this difference signifi- 
cant? If we restrict the fit to low redshift [z < 0.3) clusters 



only, we find a best-fit value of K 



370 



+ 130 



keV cm 2 



{\l = 9.69/12 = 0.81). This is consistent with the best-fit 
for the whole sample. Fitting only high redshift (z > 0.3) 
clusters we find a best-fit value of K — 220^o° keV cm 2 
(xl = 7.64/5 = 1.53). Therefore, there is only a marginal 
statistical difference between the best-fit values of K from 
the low and high redshift clusters. A very similar result is 
deduced from an examination of the S v ^ a rc/ fv — Tx rela- 
tion (which is not as sensitive as the yo~Tx relation to the 
entropy floor; MBHB03), where the low redshift clusters 
are best fit by entropy floor models with K ~ 300 keV 
cm 2 while the high redshift clusters are best fit by entropy 
floor models with K <~ 200 keV cm 2 . 



3.5. The SZ effect— Lx Relations 

The theoretical yo — Lx relations have the interesting 
properties that they do not evolve strongly with redshift 
and, more importantly, are virtually insensitive to the 
presence of an entropy floor. The reason why this hap- 
pens to be the case is not because the individual quan- 
tities, yo and Lx, are unaffected by an entropy floor (on 
the contrary, they are greatly modified), but because they 
are both affected in a very similar manner. Therefore, 
this correlation is less than ideal when it comes to probing 
the non-gravitational entropy of galaxy clusters. However, 
the relation is still of great interest because it provides a 
valuable consistency check of the other scaling relations 
studied here. 

Figure 6 is a plot of the observed yo — Lx relation. This 
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Fig. 7. Residual plots for the S v , a rcl ' fv — Lx relation. Left: 
Residuals from a comparison between the entropy floor model with 
Ko = 310 keV cm 2 and the observational data (x 2 = 4.27). Neglect- 
ing the two massive "cooling flow" clusters, A1835 and RXJ1347.5- 
1145, we converge on Ko = 385 keV cm 2 as the best-fit model with 
X 2 = 1.15. Right: Residuals from a comparison between the self- 
similar model and the observational data (x 2 = 17.30). The symbols 
have the same meaning as in Figure 1. For clarity, error bars for the 
abscissa arc not displayed. 

is superimposed on the predicted z = 0.2 (thick lines) and 
z = 0.5 (thin lines) relations for the self-similar model 
(dotted) and the K = 100 (short-dashed), 300 (long- 
dashed), 500 (dot-dashed), and 700 keV cm 2 (solid) en- 
tropy floor models. Since the theoretical yo — Lx relations 
do not evolve strongly with redshift (compare the thin and 
thick lines), a qualitative "by eye" comparison is possible. 

A clear correlation between a cluster's central Compton 
parameter and its bolometric X-ray luminosity is appar- 
ent in Figure 6. The relation is quite tight and follows the 
predicted trends. Even the two massive cooling flow clus- 
ters A1835 and RXJ1347.5-1145 seem to follow the pre- 
dicted trends. Thus, the observed yo — Lx relation gives 
us confidence that our general understanding of ICM is 
basically correct. That being said, the quality of the fit 
is not great in a statistical sense (xl = 67.36/16 = 4.21). 
The high value of xl undoubtedly arises from the random 
scatter present in the observed relation. The presence of 
this scatter may be inconsistent with theoretical predic- 
tions. We note, however, that observational systematic 
uncertainties for yo and Lx are on the order of 10% each 
(Reese et al. 2002). In addition, Mushotzky & Scharf 
(1997) have found that estimates of the X-ray luminosities 
of clusters can vary from study to study by up to 20%, 
most likely attributable to different measurement tech- 
niques. Conservatively estimating the total uncertainty 
on Lx to be 20%, we find a significantly improved fit of 
xl = 19.76/16 = 1.24. As expected, the entropy floor level 
is not constrained by the data. 

While the y — Lx relation is insensitive to the en- 
tropy floor level, the S v . arc / f v — Lx relation is not. Fit- 
ting all 15 clusters for which we have estimates of both 
the luminosity and SZ effect flux density, we find a best- 
fit entropy floor level of K = 310^70 keV cm 2 . How- 
ever, the fit is not a good one, as is evident from the 
residuals plotted in Figure 7 and the calculated reduced- 
X 2 (xl = 34.34/14 = 4.27). The two massive cooling 
flow clusters A1835 and RXJ1347.5-1145 arc obvious out- 
liers. Ignoring these two clusters, we obtain K = 385^^ 
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keV cm 2 (2a error bars) and a significantly improved fit 
(xl = 13.79/12 = 1.15). This best-fit value of K is con- 
sistent with the results of §3.2 and §3.3 and also with X-ray 
observations of nearby massive clusters. Splitting the sam- 
ple into two redshift bins (< 0.3 and > 0.3), we also find 
there to be no difference in the entropy floors of "nearby" 
and "distant" galaxy clusters. This is the same as was 
found for the S v ^ arc j f v — yo relation. 



3.6. Summary of Scaling Relations 

Every single SZ effect scaling relation that we have ex- 
amined is consistent with or requires a high value for the 
entropy floor level, Ko- In fact, several of the trends, such 
as the y - T x , yo - M(r 500 ), and S u . arc / f v - y relations, 
rule out the standard self-similar model at many sigma. 
Neither of the relations show any convincing evidence for 
strong evolution in K out to the limit to which our sample 
extends (z <~ 0.7). 

It is interesting that the estimates of K from the vari- 
ous relations do not always agree. For example, the best- 
fit entropy floors from the yo — Tx and S u , arc / /„ — Tx 
trends are consistent with results from studies of X-ray 
scaling relations of nearby massive clusters (e.g., Babul et 
al. 2002) and the results of our SZ effect-luminosity rela- 
tions but are marginally lower than the results from our 
Su.arc/ fv — yo and yo — M(r§oo) relations. A conservative 
estimate of the true value of K , however, must fall in be- 
tween the results of each of the individual relations; i.e., 
300 keV cm 2 < K < 600 keV cm 2 . This is clearly illus- 
trated in Figure 8 through a plot of A% 2 vs. Ko for the 
relations we have examined (A\ 2 = x 2 — Xbf where Xbf 
is the x 2 f° r the best-fit power-law model for each of the 
scaling relations). A naive simultaneous fit to all of the 
relations, ignoring correlated variables and errors, yields 
Ko = 430rJ" keV cm 2 (99% level), which is remarkably 
similar to that derived solely from X-ray data (c.f. Babul 
et al. 2002). 

At present, it is unclear why some of the scaling relations 
do not converge on the same value of Ko- A more detailed 
analysis of this would require taking into account all of the 
observational and analysis biases (associated with both X- 
ray and SZ effect), which is beyond the scope of this paper. 
We expect that it will be possible to get a much more firm 
handle on these differences in the near future. First, the 
number of high redshift clusters observed through the SZ 
effect is increasingly rapidly (see Reese et al. 2002) and 
with it comes improved statistics. Second, and perhaps 
more important, a number of new experiments (or sub- 
stantial upgrades to existing ones) are already under con- 
struction and will greatly improve the quality of the SZ 
effect observations. For example, the bandwidths of the 
future SZA and the upgraded OVRO array are expected 
to be nearly an order of magnitude larger than that of 
current interferometers and, as such, will lead to substan- 
tial improvements in the signal-to-noise ratios of SZ effect 
data. Tighter constraints on the SZ effect surface bright- 
ness profiles of clusters will then be possible and, in turn, 
more stringent limits can be placed on the entropy distri- 
butions of these clusters. Indeed, we demonstrate below 
that observations with the SZA and the upgraded OVRO 
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Fig. 8. Constraints on Ko ■ The dot-short-dashed, dot-long-dashed, 
dashed, and dotted lines represent the yo — Tx, Si/, arc/ fv — Lx, 
yo — M(r5oo), and S Vl arc/fv — yo relations (the four most entropy 
floor-sensitive relations), respectively. The cooling flow clusters 
A1835 and RXJ1347.5-1145 have been excluded from the fits. 

array will allow one to constrain the amount of excess en- 
tropy in clusters all the way out to z ~ 2 and, best of all, 
without the need for any X-ray results. 

4. THE FUTURE: OBSERVATIONS WITH SZA/OVRO 

It was briefly discussed in §3.1 that current SZ effect 
data alone cannot tightly constrain the three parameters 
of the surface brightness (3 model. A large degeneracy 
is present in the shape parameters, (3 and 9 C (e.g., Carl- 
strom et al. 1996; Grego et al. 2000; 2001). A number 
of authors have dealt with this problem by using the val- 
ues of (3 and 9 C obtained from modeling the X-ray surface 
brightness of the cluster (e.g., Patel et al. 2000; Graingc 
et al. 2002a; 2002b; LaRoque et al. 2003). In the case of 
low/intermediate redshift clusters, X-ray observations pro- 
vide tighter constraints on these parameters than do SZ 
effect observations. However, a truly X-ray-indepcndcnt 
probe of the ICM requires that only SZ effect data be 
used in the analysis. By using a joint maximum-likelihood 
fit to both the X-ray and SZ effect surface brightnesses of 
clusters, Reese et al. (2000; 2002) found a compromise be- 
tween these two scenarios. Pushing the analysis of the ICM 
to higher redshifts than considered by Reese et al. (2000; 
2002), however, will likely bring us into a regime that is 
uniquely accessible to SZ effect observations. Fortunately, 
substantial upgrades are planned for currently operational 
arrays [for example, the current BIMA and OVRO arrays 
are being upgraded and merged into the Combined Array 
for Research in Millimeter- wave Astronomy (CARMA) 6 ] 
and a number of new interferometers are planned as well 
(e.g., SZA 7 , AMiBA, and AMI; Holder et al. 2000; Lo 
et al. 2000; Kncissl et al. 2001). These experiments will 
substantially improve the quality of SZ effect observations. 
Below, we compare the expected performance of the SZA 
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and the (soon to be) upgraded OVRO array with the cur- 
rently operational BIMA/OVRO arrays. We show that 
the degeneracy in the — 9 C plane is greatly reduced with 
data from the SZA or the upgraded OVRO array and this 
will allow for accurate determinations of yo and S v and 
without the need for any X-ray results. 

4.1. Mock Observations 

To compare the arrays, we first generate "mock" ob- 
servations of model galaxy clusters. We follow a method 
similar to that outlined by Holder et al. (2000). First, 
we create a Compton parameter map, y(9 x ,6 y ), for each 
model cluster. These maps are converted into SZ effect 
intensity decrement maps. Before we can do this, how- 
ever, we must assume an observing frequency for each of 
the interferometers. The current BIMA and OVRO arrays 
use amplifiers which operate over the range 26 — 36 GHz 
(Reese et al. 2000). We assume a frequency centered on 
30 GHz for these arrays. The upgraded OVRO array and 
the SZA are expected to have amplifiers which operate at 
26 — 36 GHz and at 85—115 GHz. We assume frequencies 
centered on 30 and 90 GHz for these arrays. Thus, we gen- 
erate six different decrement maps for each model cluster: 
30 GHz maps for the current BIMA/OVRO arrays and 30 
and 90 GHz maps for the SZA and upgraded OVRO array. 

Interferometers do not image the sky. Rather, they mea- 
sure the Fourier transform of the SZ effect intensity decre- 
ment multiplied by the primary beam. A particular value 
of this observable is referred to as a "visibility" and is given 

by 

v(u,v)=j v [ y (e x ,e y )A(o x ,e y )e 2 ^ ue *+ ve *Ue x de y (8) 

where u and v are the conjugate variables, 9 X and 9 V are 
the projected sky coordinates, and A(0 X , 9 y ) is the primary 
beam. At 30 GHz, the sensitivity patterns for the primary 
beams of the BIMA and OVRO arrays (both current and 
upgraded) and the SZA are, or will be, well-represented 
by Gaussians with FWHM ss 6.6', 4.2', and 10 - 12', re- 
spectively. We assume a FWHM = 10.8' for the SZA. For 
more on the above formalism see White et al. (1999). 

We multiply the SZ effect maps by the Gaussian primary 
beams described above and Fourier transform t he result . 
We discard all visibilities for baselines (R UtV = V« 2 + v 2 ) 
that are not probed by these interferometers. At 30 GHz, 
the BIMA and OVRO arrays probe multipole moments 
of £ = 2nR u>v > 4000 and 7000, respectively. We as- 
sume uniform coverage over the ranges 4200 < I < 20000 
and 7200 < I < 20000 for these arrays, respectively. The 
smaller 3.5 m dishes of the SZA will be able to probe 
shorter baselines (larger angular scales) and are expected 
to sample £ > 2000. We assume uniform coverage over the 
range 2400 < i < 20000 (also at 30 GHz). 

Finally, we add random Gaussian noise to each of the 
mock observations. At 30 GHz, we assume system temper- 
atures (scaled to the atmosphere) of 45 K for the BIMA 
and OVRO arrays and 30 K for the SZA. At 90 GHz, we 
assume system temperatures of 100 K for the SZA and for 
the upgraded OVRO array. The current BIMA/OVRO 



arrays have effective noise bandwidths of approximately 
550/1000 MHz, while both the upgraded OVRO array and 
the SZA are expected to have bandwidths of 8 GHz. A to- 
tal integration time of 40 hours is assumed for each of the 
arrays at both 30 and 90 GHz. For specificity, the mock 
observations for the current BIMA and OVRO arrays have 
rms noise levels of approximately 150 /xJy beam -1 and 70 
/zJy beam -1 , respectively, while the upgraded OVRO ar- 
ray has noise levels of 25 /uJy beam -1 and 50 /zJy beam -1 
at 30 and 90 GHz (respectively) and the SZA has noise 
levels of 40 /xJy beam -1 and 140 /xJy beam -1 at 30 and 
90 GHz, respectively. 



4.2. Analysing the Mock Observations 

To analyse the mock observations, we follow the method 
described in a number of observational papers based on 
genuine BIMA/OVRO data (e.g., Carlstrom et al. 1996; 
Grego et al. 2000; 2001; Reese et al. 2000; 2002; Joy et 
al. 2001; LaRoque et al. 2003). We model the observa- 
tions with the isothermal models 8 . We create model SZ 
effect maps for for various choices of the three free param- 
eters (yo, 9 C , and 0), multiply the maps by the appropriate 
primary beams, Fourier transform the maps, and compare 
the results to the mock observations in §4.1 via a x 2 statis- 
tic. The best-fit parameters are those which result in the 
minimum value of x 2 • Note that the comparisons between 
the models and the mock observations are done in the 
Fourier domain (also referred to as the u — v plane) , where 
the noise characteristics and spatial filtering of the inter- 
ferometers arc well understood. Like the observers, we do 
not "deconvolve" the mock observations for analysis. 

To demonstrate the quality of the data to be produced 
by the SZA and the upgraded OVRO array, we plot in Fig- 
ure 9 confidence contours for the and 9 C parameters from 
"observations" of a T x = 6.7 keV (M tot w 5.6 x 1O 14 M ) 
cluster at z = 1. The third free parameter, yo, is allowed 



(i.e., 
onto 

AX 2 
The 
68% 
The 



to assume its best-fit value at each pair of (3 and 9 C 
the plot is a projection of the confidence "volume" 
the — 9 C plane). The contours correspond to a 
and the filled squares indicate the best-fit models, 
projection of A% 2 = 1 line onto the axes gives the 
confidence interval for each of the two parameters, 
contour labeled "current" is the result of fitting the SZ 
effect surface brightness model to the mock observations 
with the currently operational BIMA/OVRO arrays (both 
arrays have similar sensitivity patterns). The contours 
labeled "SZA" and "OVRO" are the results of simultane- 
ously fitting the surface brightness model to the mock 30 
and 90 GHz SZA observations and the mock 30 and 90 
GHz upgraded OVRO observations, respectively. Lastly, 
the contour labeled "SZA+OVRO" is the result of simul- 
taneously fitting all four upgraded OVRO and SZA mock 
observations. 

First, it is obvious that the "current" mock observa- 
tional data do not tightly constrain and 9 C and, further- 
more, these parameters are strongly correlated, as found 
when modeling genuine BIMA/OVRO data (see, e.g., Fig. 
3 of Grego et al. 2000; Fig. 2 of Grego et al. 2001). This 



6 For information on the CARMA see http://www.mmarray.org 

7 For information on the SZA sec http://astro.uchicago.edu/sze/survey.html 

8 Note. — we do not assume that the model cluster used to make the mock observations is necessarily isothermal. 
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Fig. 9. Comparison of constraints on the surface brightness profile 
of a distant, massive cluster. The filled squares indicate the best-fit 
parameters from modeling "mock" observations of the model cluster 
(the square at /3 = 1.74 is for "current" data while the square at 
fi = 1.59 is for "SZA" data). The contours give the 68% confidence 
regions on the two parameters (Ax 2 = 1). The normalization, yo, 
is allowed to assume its best-fit value at each point in the plot (see 
text). 

gives us confidence that we have produced realistic mock 
observations. The large degeneracy in the j3 — 8 C plane 
is why a number authors elect to use X-ray constraints 
on these parameters instead or, in the case of Reese et 
al. (2000; 2002), use a simultaneous fit to both X-ray 
and SZ effect data. However, data from the upcoming 
SZA and the upgraded OVRO array will be able to place 
much tighter constraints on these parameters (compare 
the "OVRO" and "SZA" contours with the "current" con- 
tour). Note that the correlation between (3 and 8 C remains 
for the SZA and OVRO "data", but its size has been dra- 
matically reduced. A simultaneous fit to both the SZA 
and upgraded OVRO data yields even better constraints 
on these parameters (shaded region). 

To get an idea of how well the future SZ effect data can 
constrain Kq, the entropy floor level, we use the 68% con- 
fidence volume (/?, C , yo) for the "SZA+OVRO" contour 
to "measure" yo and S v ^ arc j '/„. The inferred statistical un- 
certainty associated with the central Compton parameter 
and integrated SZ effect flux density within the central 1' 
for this cluster is only about 10% and 15%, respectively. 
Comparing this to the predicted S v . arc / j v — yo relations 
(Figure 2 of MBHB03), it should, therefore, be possible to 
constrain the entropy floor level of this cluster to within 
50-75 keV cm 2 or so. This is comparable with statisti- 
cal uncertainties associated with X-ray measurements of 
nearby clusters. This is remarkable considering that no 
X-ray "data" was used in the analysis and the cluster lies 
at z = 1. We also find that reasonably accurate measure- 
ments of Ko are possible for clusters all the way out to 
z ~ 2. 



The SZA and upgraded OVRO array will be excellent 
tools for probing the non-gravitational entropy of distant 
clusters. Because X-ray data will not be required to con- 
strain the shapes of the SZ effect surface brightness pro- 
files of clusters observed with these planned interferom- 
eters, comparisons of data to predicted scalings (such as 
S v ,arc/ fv — yo) will provide independent constraints on the 
properties of the intracluster gas. It will then also be pos- 
sible to take advantage of the redshift-independence of the 
SZ effect and monitor the evolution of non-gravitational 
processes in clusters right back to the epoch of cluster for- 
mation itself. 

5. DISCUSSION 

Up until now, measurements (both direct and indirect) 
of the entropy floors of massive galaxy clusters have been 
limited to X-ray observations. Furthermore, these past 
X-ray studies have generally focused on nearby [z <~ 0) 
clusters and, as such, little is known about the evolu- 
tion of the entropy floor (and the non-gravitational pro- 
cesses that produce it) with cosmic time. In the compan- 
ion paper, we explored the extent to which the thermal 
Sunyacv-Zcldovich effect is modified by the presence of 
an entropy floor. Because it depends differently on the 
temperature and density of ICM and, also, because it is 
redshift-independent, the SZ effect could potentially be a 
very powerful, independent test of the entropy floors of 
even the most distant galaxy clusters. 

The central focus of the present paper was to com- 
pare our theoretical relations from MBHB03 (including 
one that can potentially be measured through SZ effect 
observations only) to available high redshift SZ effect data 
from the literature to determine if the SZ effect data sup- 
port the presence of an entropy floor and, if so, how does 
the inferred level of that floor compare with that required 
to match local X-ray trends. This is the first time such a 
comparison has been done and we have made use of the 
largest compilation of high z SZ effect clusters to date. 
A detailed analysis of seven different SZ effect scaling re- 
lations indicates that the entropy floor in clusters with 
0-14 < z < 0.78 is between 300 keV cm 2 < K < 600 keV 
cm 2 and there are no strong indications for evolution in 
Ko over that redshift interval. Our estimate for the value 
of Ko is remarkably similar to that derived from studies 
of X-ray scaling relations of nearby (z ~ 0) massive clus- 
ters, which suggest that K > 300 keV cm 2 (e.g., Tozzi & 
Norman 2001; Babul et al. 2002: McCarthy et al. 2002). 

At present, the source (or sources) of the "excess" en- 
tropy is still not known. What constraints can be placed 
on the possible sources by the results of the present study? 
First, as in previous studies of X-ray scaling relations, 
our analysis indicates that the entropy of the ICM has 
been significantly raised by some non-gravitational pro- 
cesses). In terms of thermal energy, this corresponds to 
a few keV per particle for massive clusters. This means 
that supernovae explosions probably cannot be the sole 
contributor to the entropy floor, since they are expected 
to impart < 0.3 keV per particle (e.g., Valageas & Silk 
1999; Balogh et al. 1999; Wu et al. 2000). This was 
previously known but the present study, which offers an 
independent examination of the ICM, reinforces this con- 
clusion. Recently, it has been speculated that quasar out- 
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flows could be source of the excess entropy (e.g., Nath & 
Roychowdhury 2002). The entropy requirements deduced 
in the present study (and previous X-ray studies) are prob- 
ably met by quasar outflows but it isn't yet known what 
mechanism (if any) couples the outflows to the ambient 
ICM. Alternatively, and somewhat paradoxically, radia- 
tive cooling has also been shown to raise the mean entropy 
of ICM. It is possible that cooling in combination with 
supernovae and/or quasar outflows could be responsible 
for the observed SZ effect and SZ effect- X-ray relations. 
Whatever the source may be, it must reproduce the fact 
that K does not change significantly out to z <~ 0.7. It 
could well be that this trend will become a critical piece 
of information for discriminating between the various the- 
oretical models currently being proposed. An unchanging 
value of Kq with redshift is obviously consistent with the 
generic "preheating" scenario; however, it remains to be 
seen whether it is consistent with more realistic heating 
models that distribute entropy non-uniformly and over an 
extended periods of time. Without a detailed analysis, 
it is difficult to say whether or not it is consistent with 
radiative cooling contributing significantly to the excess 
entropy. We are currently in the process of examining the 
effects of radiative cooling on SZ effect scaling relations 
(c.f. the discussion in the companion paper, MBHB03). 

Current SZ effect data cannot tightly constrain the sur- 
face brightness profiles of clusters. This prevents the SZ 
effect from being used as an independent (of X-ray) probe 
of the entropy floor and the ICM in general. Thus, any 
advantages that the SZ effect has over the X-ray emission 
(e.g., redshift independence) are severely diminished be- 
cause X-ray data is needed to help constrain the shape 
of the surface brightness profiles. Therefore, an additional 
aim was to examine the ability of the next generation of SZ 
effect experiments to probe non-gravitational entropy in 
distant clusters. We have shown that the SZA and the up- 
graded OVRO array will produce high signal-to-noise data 
that will allow one to tightly constrain the surface bright- 
ness profiles of even very distant clusters and without the 
need for any X-ray results. As in the present study, these 



surface brightness profiles can then be compared to theo- 
retical predictions in order to place stringent constraints 
on the level of the entropy floor. It will be very interesting 
to see if the trend of constant K (with redshift) deduced 
here holds up and, if so, to determine how far back in 
redshift it extends. 

Aside from the appearance and structure of individual 
clusters, our work has implications for universal SZ effect 
quantities, such as the SZ effect angular power spectrum, 
SZ effect cluster source counts, and the mean Compton 
parameter of the universe. These quantities can be used 
to measure cosmological parameters and test cluster for- 
mation scenarios (see Carlstrom, Holder, & Reese 2002 
for a comprehensive review). A number of studies have 
already examined how non-gravitational gas physics mod- 
ifies these quantities (e.g., Holder & Carlstrom 2001; Cav- 
aliere & Menci 2001; da Silva ct al. 2001; Springel et al. 
2001; White et al. 2002). However, none of these studies 
have invoked entropy injection at the high level estimated 
in the present analysis. Generally, low levels of entropy 
injection (K ~ 100 keV cm 2 ), which are consistent with 
X-ray measurements from low mass groups (Ponman et 
al. 1999), were implemented. We are in the process of 
investigating how the power spectrum, source counts, and 
mean Compton parameter are modified by higher initial 
entropies. 
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sz effect observational data 
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Note. — S v . arc l '/„ expressed in mJy. 
a Reese et al. (2002) 
b Joncs et al. (2003) 
c Holzapfel (1996) 
d Holzapfel et al. (1997) 
°Lamarre et al. (1998) 
'Saunders et al. (1999) 
s Mauskopf et al. (2000) 
h Pointecouteau et al. (2002) 
'Pointecouteau et al. (2001) 
JKomatsu et al. (1999) 
k Grainge et al. (2002a) 
'Hughes & Birkinshaw (1998) 
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